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Abstract 

We study the quantum mechanics of the derivative nonlinear Schrodinger 
equation which has appeared in many areas of physics and is known to be 
classically integrable. We find that the iV-body quantum problem is ex- 
actly solvable with both bound states (with an upper bound on the particle 
number) and scattering states. Quantization provides an alternative way to 
understand various features of the classical model, such as chiral solitons and 
two-soliton scattering. 
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The derivative nonlinear Schrodinger equation (DNLS) in one dimension 
has historically found applications in many areas of physics, one example 
being circularly polarized nonlinear Alfven waves in a plasma |IJ. Related 
models have recently received fresh attention in the context of chiral Lut- 
tinger liquids; some of these models can be obtained by a dimensional reduc- 
tion of a Chern-Simons model defined in two dimensions [2-5]. The DNLS 
has some peculiarities, for instance, it is not Galilean invariant, and it has 
classical solitons which have an upper bound on the particle number and are 
chiral (with a particular sign of the momentum). In these respects, the DNLS 
differs from the usual nonlinear Schrodinger equation (NLS), although both 
of them are classically integrable (see Refs. |], 0] and references therein). 
Unlike the DNLS, the usual NLS has been studied semiclassically and quan- 
tum mechanically in great detail || [|; being equivalent to a Bose gas with 
a two-body attractive 5-function interaction |T0, ITT] . These studies have led 



to an alternative understanding of various features of the classical NLS, such 
as the solitons and their scattering. 

It therefore seems interesting to similarly analyze a quantum mechanical 
version of the DNLS. We do so in this Letter generalizing the analysis of 
the two-body problem given in Refs. [2-5]. Besides providing us with a new 
exactly solvable iV-body quantum mechanical problem, our study leads to 
a different and perhaps simpler way of deriving various properties of the 
DNLS, such as the peculiar features of the solitons and the time delay in the 
scattering of two solitons. 

We begin by considering a general Lagrangian density of the form 

h 2 

£ = ~ 77- [ d x ij>*d x 1> + i\p(ip*d x ip - d x il>*ij>) + 2pp 2 + vp z ] , (1) 

2m 

where p = ip*ip is the density with dimensions of inverse length. We have 
introduced Planck's constant h in Eq. ([!]) for later convenience, although we 
will first discuss classical mechanics. The usual NLS is obtained by putting 
A = v — in (Q), while p — v — produces the DNLS. Note that the A term 
is not invariant under Galilean transformations [|3], |J , and it flips sign under 
parity (x — > —x). It is therefore sufficient to consider the case A > 0. The 
parameters A and v are dimensionless while p has the dimensions of inverse 
length; the system is scale invariant if we set p = 0. 

We will only examine the case v — for the following reasons. Firstly, 
it is known that the Euler-Lagrange equations of motion which follow from 
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Eq. ([[]) are classically integrable if u = 0, regardless of the values of A 
and /i |, |. Secondly, it is not clear to us how to handle the up 3 term 
in (JJ) quantum mechanically. This piece would lead to a term like p 2 ^> in 
the Schrodinger equation, where \1/ is the wave function. If the positions of 
the particles are denoted by Xi,x 2 , . . . ,x^, the density operator is given by 
p{xi) = 2~2j S(xi—Xj). Hence p 2 contains highly singular terms like <5 2 (xj — Xj) 
as well as three-body terms like 5{xi — Xj)5(xi — Xk)- (A possible way of 
interpreting o ( the limit of a two-body (^-function interaction 

with infinite strength; this forces the wave function to vanish whenever the 
coordinates of two particles coincide || . However the problem then becomes 
that of hard core bosons or free fermions in one dimension which is easily 
solvable). We can avoid such singular interactions by setting v = 0. This 
gives us the DNLS with a slight generalization if p is nonzero. 

For simplicity, we will set p = in Eq. ([TJ) to start with. Later we will 
reintroduce p and will comment on an interesting phenomenon which occurs 
in that case. We now briefly summarize the classical one-soliton solution of 
([I]). The equation of motion is 

in w = - Jf- flgv + — p d x * . (2) 

2m m 

The density satisfies the continuity equation d t p + d x j = 0, where 

it) f> \ 

J = - 7T - - —P 2 ■ ( 3 ) 

2m m 

Since the system is integrable, there are an infinite number of conserved 
quantities [B] . Three of these are the particle number, momentum and energy, 



N = I dx p 

ih 

~2 



p = - l — I dx - d x r^P) 



h 2 r 

E = — J dx[ d x iP*d x il> + iXp(<ip*d x ip ~ d x r^) } ■ (4) 
With A > 0, the one-soliton solution takes the form || [| 

ip(x, t) = J p(x,t) exp [ i (kx - ut + A / dy p(y, t) ) ] , 

J —oo 
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a 2 

P ^'^ ~ 2A( V« 2 + 4F cosha(x-ft) + 2k) ' 

mt) mv 2 ha 2 

k = — — and = — . (5) 

% 2h 8m v ' 

These expressions contain two independent parameters, the velocity v and the 
inverse width of the soliton a. On computing the three conserved quantities, 
we find that < iVA < ir, and 



cosiVA 



2k 



Va 2 + 4k 2 



M = ^ taniVA , (6) 
A 

with P = Mv and E = P 2 /2M. We observe that P is always positive, 
although M, v and E are positive if < iVA < ti/2 and negative if n/2 < 
cos NX < n. (P would have been negative if we had chosen A < 0). 

We now study this system quantum mechanically, and show that the 
classical picture is recovered in the limit N — > oo. The Schrodinger equation 
describing N identical bosons is easily derived 0, |^| after defining the wave 
function 

*( Xl ,x 2 ,...,x N ,t) = (0\if(x 1 ,t)if(x2,t)...if(x N ,t)\N) , (7) 

where $(x, t) denotes the second-quantized bosonic field operator. We then 
obtain the equation 

ih d& = - f - E 91* + — E 8( Xi - Xj )(d Xi + d Xj )V , (8) 
2m m ~r 

where we have normal ordered the Hamiltonian to eliminate divergent self- 
interaction terms like 5(xi — Xj) = 5(0). We thus find a two-body 5-function 
interaction whose strength depends on the total momentum of the two par- 
ticles. Note that we do not need to worry about the relative ordering of 
d Xl + 9 X . and 5(xi — Xj) since they commute; thus the 5-function interaction 
does not affect the total momentum of the two particles. 

Before proceeding further, it is instructive to examine the two-particle 
problem. For the configuration x\ < x% y let us assume a stationary solution 
of the form 

% % 
\1> = exp [ ^( c i x i + C2X2) ] + a exp [ ^( c 2^i + cix 2 ) ] , (9) 
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where c\ and c 2 may be complex, and a will be called the reflection amplitude. 
The wave function for x 2 < %i is then obtained by Bose symmetry. For 
convenience, we will refer to the exponentials in Eq. as waves and the c n 's 
as the particle momenta even if the C n S £1X6 not real. The total momentum 
of (H) is P — c\ + c 2 . We can now go to the center of mass and relative 
coordinates X = (x 1 + x 2 )/2 and x = x\ — x 2 , and factor out the center of 
mass wave function. The relative motion is described by the equation 

h 2 h 2 \P P 2 

ihdtV = [ d 2 x - 5{x) + — ] * . (10) 

We then find that the amplitude a in Eq. (§) is given by 

c 2 - ci - i\(c 2 + ci) 



C 2 - Cl + «A(C 2 + Cl; 



(11) 



Clearly, there is a bound state if a = and P > 0. 

We now seek a bound state solution of the iV-particle Schrodinger equa- 
tion(§, 

* = eX P [ i C ^ X n } (12) 

for the configuration x\ < x 2 < . . . < xn] the wave function for all other 



configurations can be deduced by Bose symmetry. We find that ( |T2| ) satis- 
fies the boundary conditions of the 5-function interactions if the reflection 
amplitude given in Eq. fllTf ) is zero for each neighboring pair of particles, 
i.e., 

c„+i - Cn = iX (c n+ i + c n ) (13) 

for n — 1, 2, . . . , N — 1. (It is because of the absence of the reflected waves 
that our wave function ( |I~2"D has only one wave, instead of iV! as in the general 
Bethe ansatz). We now demand that the total momentum should be equal 
to P = J2n c n, and discover that the c n 's are given by 

c n = exp [i(2n-N-l)6}, (14) 

where 

exp (i2B) = . (15) 
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The parameter 9 lies in the range [0, tt/2] for A > 0. The energy is given by 
E = J2n cl/2m = P 2 /2M, where 

tan iV6» 

M = m . (16) 

tan 9 v ' 

We now ask, when do Eqs. ([l|) and (0) describe a bound state, i.e., a 
state which is normalizable if we use translation invariance to fix the centre 
of mass at some particular place? It is easy to show that the necessary and 
sufficient condition for this is that 



sin 9 
sin 



CUI1U1L1UI1 1UI liiis is mat 

a9 ' , sm(N-l)9 sinW , . 

— g S in(JV + 1 - 2n)B = P ^—^ > , (17) 

for all values of I from 1 to N/2 if N is even and to (N — l)/2 if N is odd. 
The conditions in (\Tf\j arise on demanding that the probability vj/*\I/ should 
go to zero if we take the / particles on the right (xn+i-i, xn+2-i, ■ ■ ■ , %n) to 
oo and the / particles on the left (xi, X2, ■ ■ ■ , %i) to —00 (thereby keeping the 
centre of mass fixed). For small values of N, we can check whether or not 
Eqs. fllTf ) are satisfied; we find an intricate pattern of allowed values of 9 
and P (which sometimes has to be negative) for which a bound state exists. 
The situation simplifies in the limit N — ► 00 and 9 — > keeping N9 fixed. 
We then find a bound state if and only if N9 < 71 and P > 0. Further, in 
this limit, 9 = X so that the classical and quantum formulae (||) and ([16]) for 
the masses agree; we therefore identify the quantum bound state with the 
classical soliton. (It would be interesting to derive the soliton profile in Eq. 
(||D from the wave function in (|T2|) using a technique given in Ref. ||). We 
thus see that only a finite number of particles can be bound for a given value 
of 9. It is quite remarkable that if n/2 < N9 < n, the energy of the bound 
state can be lowered arbitrarily by giving it sufficiently large momentum. 

Next we study the scattering of bound states. As a warmup exercise, 
consider the scattering of a bound object of N — 1 particles with momentum 
P — p with one particle of momentum p, where p is real. For x\ < . . . < x^, 
we find that the wave function can be written as a superposition of N waves, 



2 

* = exp [-(cixi + . . . + C7v_i£at_i + px N )] 
n 

+Aexp [-{px\ + c\x 2 + • • • + c N -ix N )} 
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N-2 i 

+ a n exp [t{c\Xi + . . . + c n x n +px n+1 + . . . + c Ar _iX iV )]{18) 

n=l 11 

where the c n 's are given by Eq. ( |14|) with iV and P replaced by N — 1 and 
P — p respectively in that equation. The first two waves in ([18]) correspond 
respectively to configurations in which the the bound state is entirely to the 
left and entirely to the right of the free particle; A is therefore the transmis- 
sion amplitude for the particle to go through the bound object. We now use 
Eq. (|TT| ) repeatedly to relate all the a n 's and A to each other. A is found to 
be a pure phase of the form A = C/C*, where 



N_1 , r, P-n 
C = TT ( e- 19 P e i(2n-N+l)e \ (1Q) 



n=l 



We now examine the general scattering of two bound objects (solitons), 
and again discover that there is only transmission and no reflection. We 
consider one object with Ni particles and momentum Pi and another object 
with N2 particles and momentum P2; let N\ + N2 = N. We introduce N\ 
numbers satisfying Eq. ( |T^) (with N, P replaced by Ni, Pi), and N 2 
numbers satisfying (O) with N2, P%. For a given configuration x\ < 
x 2 < ... < xn, we then find that the wave function is a superposition of 
several waves, each wave having the particle momenta as some permutation 
of the numbers and We do not get all the N\ possible permutations 
due to the absence of reflection within the momenta c^'s, or within the 
c^'s, as discussed around Eq. (|T3|). The permutations allowed are those in 



which particles i and j can have the momenta and c n l_i (or momenta 

c^f 1 and Cn+i) only if i < j. The number of such permutations is N\/Ni\N 2 \- 
Two of these describe configurations in which the object with Ni particles 
is entirely to the left (or to the right) of the object with A^2 particles; all 
the other permutations describe configurations in which the Ni particles are 
interspersed amongst the A^ 2 particles. Thus the wave function has the form 

i Ni N 2 

* = exp [ -(53 C^X n + Cn )x N 1+ n) ] 

n n=l n=l 

+ A exp [ -<X c^x n + ^ c^x N2+n ) ] 

n n=l n=l 

+ all the other waves , (20) 
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where A denotes the transmission amplitude for one bound state to pass 
through the other. We compute A by successively passing each of the A2 
particles on the right through each of the iV\ particles on the left, and using 
the expressions in Eq. ([11]) at each such crossing. We finally discover that 
A = exp(i25) = C/C*, where 

N 2 / P P 

Q _ TT ( r<2 i(2n~2-N 2 )e _ f_l ANxt 

*i 2 V s inAr 2 ^ sin Ntf 

Ni p p 

. TT ( r2 e -iN20 _ t± i(2n-JV x )0 \ ^l) 

^ V s inA^ sin N x > ' V ' 

Let us consider the case in which the velocities Vi = (2E i /M i ) 1 ^ 2 satisfy 
V\ > v 2 > 0. For weak coupling {0 — > with Aq, AT 2 held fixed, so that 
Mi = mNi), we find the phase shift 



' Ei /Mi + JE 2 /M 2 
25 = 26NiN 2 \ ' V . 1 + O(0 2 ) . (22) 
I E1/M1 - JE 2 /M 2 



From this, we can compute the time delay suffered by the bound state with 
energy E x ; thus At x = 2hd5/dEi |12[. The expression in (22) and the 



resultant time delay Ati agree with the results in Ref. || if we set Ni = N 2 . 

Let us now include the term proportional to \i in Eq. (]!]). This adds 
(2h 2 n/m) J2i<j ${ x i ~ x j) to the right hand side of the Schrodinger equation 
(ffl). For a bound state of iV particles with momentum P, we find that the 
numbers c n in Eq. (|12D are given by 



Nhu\ sm9 . .. , T ha 

TT ) . 777; exp i(2n - JV - 1)0 + — . 23 

tanfl' sinAT0 H 1 v 11 tan0 v ; 



As before, we find that in the limit 0^0 with NO held fixed, the bound 
state exists only if NO < it; in addition, we need P > Nhfi/ tan (9. More 
interestingly, we observe that the energy E = J2n c n/^ m * s minimum at a 
nonzero value of the momentum given by 

tan0 V Ntan0 / v ; 

if fitemNO < 0. For weak coupling (A^0 < 7r/2), we can understand this 
result as follows. We know that there is a zero momentum soliton if \x < 
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(attractive interaction) even if 9 = [H, O]; the energy of this is given by 



h 2 a 2 

E = - -f- (N 3 -N) , (25) 
om 



as we can see by setting P = and taking the limit 9 — > in Eq. (0). Now 
if # is small and positive, the strength of the attractive 5-function interaction 
is increased if all the particles move with positive momentum; this lowers the 
energy by effectively increasing the value of /z 2 in (j25P . Thus zero momentum 
is not the state of lowest energy if 9 is nonzero and \x < 0. 

For completeness, we would like to mention the states in which all the 
particles have real momenta p n . (However these purely scattering states are 
not the lowest energy states of our system). It is convenient to put the system 
on a circle with circumference L, and consider a particular ordering of the 
positions < X\ < . . . < < L. The wave function is then given by the 
general Bethe ansatz with a superposition of iV! waves. Following Ref. [|TH|] . 
we impose periodic boundary conditions 

W(p,x 2 ,---,x N ) = V(x 2 ,...,x N ,L) . (26) 

We then find that the p n 's are related to each other by the N equations 

T N — 

2nl n + ir(N- 1) + 2 £ tan^ 1 ( Pn . Pl ) , (27) 



where the J n 's are integers. These equations may be solved numerically. 

Finally, we ask whether the system has a well-defined ground state in the 
thermodynamic limit N, L — > oo keeping N/ L fixed. Our earlier analysis 
indicates that the answer is no, even if /i is zero or even positive (repulsive). 
We have seen that a bound state with iV particles can have an arbitrarily 
low energy if its momentum is large and 7r/2 < N9 < it. UN is very large, 
the system can lower its energy arbitrarily by forming a number of large 
momentum bound objects with particle numbers N\, N 2 , . . . , (adding upto 
N), such that tt/2 < N { 9 < ir for i = 1, 2, . . . , k. 

To conclude, the quantization of the DNLS has produced a rich structure. 
It would be interesting to consider other classically integrable systems and 
see if they can be quantized in order to shed new light on them. 

I thank Raj at Bhaduri for discussions and the Department of Physics and 
Astronomy, McMaster University for its hospitality during the course of this 
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work. This research was supported by the Natural Sciences and Engineering 
Research Council of Canada. 



Note Added: 

After writing this paper, I learnt that similar work has been published 
earlier O, |I3fl; I thank A. G. Shnirman for pointing this out. 
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